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Abstract. We have shown that the cluster-mass reconstruction method which combines strong and weak grav- 
itational lensing data, developed in the first paper in the series, successfully reconstructs the mass distribu- 
• tion of a simulated cluster. In this paper we apply the method to the ground-based high-quality multi-colour 

data of RX J1347.5— 1145, the most X-ray luminous cluster to date. A new analysis of the cluster core on 
very deep, multi-colour data analysis of VLT/FORS data reveals many more arc candidates than previously 
known for this cluster. The combined strong and weak lensing reconstruction confirms that the cluster is indeed 
very massive. If the redshift and identification of the multiple-image system as well as the redshift estimates of 
the source galaxies used for weak lensing are correct, we determine the enclosed cluster mass in a cylinder to 
CLf M(< 360h _1 kpc) = (1.2±0.3) x 10 15 M o . In addition the reconstructed mass distribution follows the distribution 

found with independent methods (X-ray measurements, SZ). With higher resolution (e.g. HST imaging data) 
5_j ■ more reliable multiple imaging information can be obtained and the reconstruction can be improved to accuracies 

' greater than what is currently possible with weak and strong lensing techniques. 
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1. Introduction have an advantage in constraining the mass at much larger 

radii than strong lensing, one of main limitations for both 

Clusters of galaxies have been a focus of a very intense on- gtrong and weak lenging ig the problem of the mass _ sheet 

going research. Especially important for many cosmologi- degeneraC y (i.e. the mass profile of the cluster can only be 

cal applications is a good determination of their mass. One determined up to a constan t). In the absence of redshift 

way to obtain their masses is to use the gravitational lens- information from ind i vidua i source s and the lens, one can 

mg information, both from multiple image systems (strong break tMg degeneracy only by making assump tions about 

lensing) as well as from distortions of background sources thc underlying potential. Different assumptions, however, 

(weak lensing). Many weak and strong lensing cluster mass can lead to discrepant results on the cluster mass . In this 

reconstructions have been successfully performed in the wQrk we therefore use individua i re dshifts of background 

past (see e.g. Clowe & Schneider 2001, 2002, Gavazzi et al. SQurces tQ overcome this prob i em . As shown in Bradac 

2004 for examples of weak lensing and e.g. Kneib et al. gt al ( 20 04a), by using these and by extending the recon- 

2003, Smith et al. 2004 for a combination of weak and struction to the inner partg of the cluster we are effectively 

strong lensing). While weak lensing mass reconstructions ab j e to break ^ b jg degeneracy 
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method in which we extend the weak lensing formalism to 
the inner parts of the cluster, use redshift information of 
the background sources and combine these with the con- 
straints from multiply imaged systems. Using simulated 
data we have shown that the method is successful in recon- 
structing the mass distribution of a cluster, and yields an 
excellent agreement between the input and reconstructed 
mass also on scales within and beyond the Einstein radius. 

Encouraged by the success of our method, we apply 
it to the weak and strong lensing data for the redshift 
0.451 cluster RX J1347.5-1145 (Schindlcr et al. 1995), 
the most X-ray luminous cluster known to date. Due 
to its record holding, this cluster has been a subject of 
many studies in X-ray (Schindler et al. 1995, 1997, Allen 
et al. 2002, Ettori et al. 2004, Gitti & Schindler 2004) 
and optical (Fischer & Tyson 1997, Sahu et al. 1998, 
Cohen & Kneib 2002, Ravindranath & Ho 2002). It has 
also been detected through the Sunyaev-Zel'dovich effect 
(Pointecouteau et al. 2001, Komatsu et al. 2001, Kitayama 
et al. 2004). Yet the mass determinations based on X- 
ray properties, SZ effect, velocity dispersion measurement, 
strong and weak lensing have all yielded discrepant results 
(see Cohen & Kneib 2002 for a summary). 

For the purpose of mass reconstruction we use 
VLT/FORS data on a field of 3.8 x 3.8 arcmin 2 in U, 
B, V, R, and I bands. We also use Ks-band data from 
VLT/ISAAC to obtain more reliable photometric redshift 
estimates. The shape measurements for the weak lensing 
reconstruction is performed on two FORS bands, R and I. 
The strong lensing properties of this cluster are analysed. 
From previous data sets five arc candidates were reported 
(Schindler et al. 1995, Sahu et al. 1998); using the new 
multi-colour data we conclude that only two possibly be- 
long to the same multiple image system. Furthermore, we 
searched for additional images belonging to this system 
and identified a third possible member. Several new arc 
candidates were found as well and are presented in this 
work. Particularly interesting is a very red arc candidate 
with two components, located at a distance of 1 arcmin 
from the brightest cluster galaxy (BCG). In addition, we 
detect further elongated structures, some of them have 
been previously indicated by Lenzen et al. (2004) who de- 
veloped and use an automated arc searching routine. 

This paper is organised as follows. In Sect. 2 we de- 
scribe the observations and give a brief outline of the data 
reduction process. In Sect. 3 we describe how we search for 
multiply imaged systems. In Sect. 4 we give the results of a 
combined strong and weak lensing reconstruction and the 
cluster luminosity measurements. We conclude in Sect. 5. 

2. Observations and data reduction process 

The optical VLT data for the current project were 
obtained with ESO proposal 67.A-0427(A-C) (P.I. S. 
Schindlcr). The data was taken with FORS1 in the high- 
resolution mode (pixel scale C'09; total field of view w 
3'.2 x 3'.2) in service mode between April and September 
2001. UBVRI Bessel filters were used in sub-arcsecond see- 



ing conditions (see Table 1 for a summary of data prop- 
erties). This allows us to estimate photometric redshifts 
for all galaxies and to support our mass and light analy- 
sis by a careful separation of foreground and background 
galaxies and cluster members (see below). Our primary 
band for the weak lensing analysis (the I) was taken in 
the 1-port read-out mode. Thus we avoid potential prob- 
lems for object shape measurements due to varying noise 
properties in the central parts of the images. For the other 
4 bands, primarily used for object photometry, the 4-port 
read-out mode was used. The data in each band consist of 
at least 20 individual exposures and were obtained with a 
dither pattern of 30'.'0 in RA and DEC in order to obtain 
clean coadded images of highest quality. 

The data reduction was carried out with a pre-release 
version of THELI, a pipeline developed specifically for 
the processing of optical single- and multi-chip cameras 
(see Schirmer et al. 2003, Erben et al. 2005); here we 
only outline our astromctric calibration which is essen- 
tial for weak lensing studies. First, we match object po- 
sitions from I-band data with those from the USNO-A2 
astrometric catalogue (Monet et al. 1998), which fixes the 
position of the individual exposures with respect to abso- 
lute sky coordinates and thus corresponds to a zero-order 
astrometric solution ("shift only"). Next, we used Mario 
Radovich's Astrometrix (see McCracken et al. 2003 and 
http : //www . na . astro . it/~radovich/WIFIX/) to fit im- 
age distortions by a two-dimensional, third-order polyno- 
mial. Hereby, the distances of the objects with coordi- 
nates in USNO-A2 catalogue and of the overlap sources 
in different images are minimised simultaneously in the 
X sense. We end up with rms residuals of w 0'.'25 for the 
USNO-A2 standard sources and « O'.'Ol for the overlap ob- 
jects. Afterwards, we extract high S/N objects from the 
coadded I-band image which are used as astrometric stan- 
dard sources (instead of USNO-A2) for the other bands. 
In all bands we achieve formally an internal astrometric 
accuracy of w O'.'Ol — 0'.'015 for the overlap sources. Most 
of the observations were done during photometric nights. 
Photometric zeropoints were deduced from the images of 
standard stars obtained as part of the standard calibra- 
tion plan of the FORS1 instrument and reduced in the 
same way as the science data. The obtained zeropoints 
are in good agreement with the general trend analysis of 
the FORS1 zeropoints. From non-photometric nights we 
only include images with a maximum absorption of 0.1 
mag in the coaddition process which is performed with 
drizzle (Fruchter & Hook 2002). 

In addition, we retrieved Ks VLT-ISAAC data (pixel 
scale C/1484; field of view « 2'.5 x 2'.5) from the ESO sci- 
ence archive (proposal ID 67.A-0095(B)). The data was 
processed with the eclipse package (see Devillard 1997). 

We create the catalogue of objects using SExtractor 
(Bcrtin & Arnouts 1996) in dual-image mode. The I-band 
image is used for detections and the images of the other 
bands are only used to measure the corresponding mag- 
nitudes. An object is considered detected if five adjacent 
pixels had a flux that exceeded the local sky noise level 
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Table 1. Properties of the data used in this work. The 5<r 
limiting magnitudes were determined with SExtractor using 
an aperture of 2". 



Filter 


Exposure time (s) 


Seeing (") 


Limiting mag. 


U 


11310 


0.97 


25.4 


B 


4800 


0.67 


26.9 


V 


4500 


0.62 


26.5 


R 


6000 


0.67 


26.3 


I 


6750 


0.57 


25.6 


Ks 


« 7200 


0.73 


21.4 



by a factor of three. All magnitudes quoted in this paper 
are in the Vega system. The photometric redshifts (using 
isophotal magnitudes for cluster members and aperture 
magnitudes for background sources; see below) of the ob- 
jects were obtained using the HyperZ package (Bolzonclla 
et al. 2000). 

2.1. Cluster member catalog 




1 2 



m y — rxij 

Fig. 1. The niB — my vs. my — mi colours for the galaxies in 
our field. Cluster members are selected to lie inside the poly- 
gon. BCG colours are given as a triangle. In addition we plot 
as crosses all the galaxies which have a photometric redshift 
estimate 0.4 < z p hot < 0.5. 

For cluster members, due to their brightness and dif- 
ferent sizes, we argue that it is best to use isophotal mag- 
nitudes to obtain accurate colour estimates. In Fig. 1 we 
plot the bib — my vs. my — mi colours for the galaxies 
in our field. To determine the colour cuts for the cluster 
member selection we first inspect the galaxies having I- 
band magnitude up to three magnitudes fainter than the 



BCG and with a distance to the BCG smaller than V . 
These are preferentially the cluster members and form a 
group around my — m\ ~ 2 and m^ — my ~ 1.8. Using this 
information we determine the following selection criteria 
for the cluster members 

0.7 (tob — my + 1) < my — m\ < 0.7 (tob — my + 1.9) 

1 < m B - my < 2 , (1) 

and we also cut out all the objects having magnitudes 
brighter than the BCG. In Fig. 1 we plot colours for all 
galaxies in our field, the BCG colours (slightly bluer than 
most other members) are given as a triangle and the poly- 
gon indicates the selection criteria we use. In addition, 
to avoid biases toward red cluster members, we add to 
the catalogue the blue galaxies with photometric redshifts 
0.4 < Zphot < 0.5 (denoted as crosses in Fig. 1). The final 
photometric redshift distribution of the cluster members 
is given in Fig. 2. The completeness of our catalogue is 
discussed in Sect. 4. 4. 

An alternative approach would be to select the cluster 
members purely from the redshift information. However 
due to uncertainties in redshift estimation the redshift dis- 
tribution of the members is relatively broad. With broad 
cuts in redshift space one can get, on the one hand, a con- 
tamination of blue, non-cluster members and on the other 
hand, some red cluster members might be missed due to 
an incorrect redshift determination. In the previously de- 
scribed method the situation is reversed. We have tested 
both selection criteria to calculate the cluster luminosities 
(see Sect. 4.4) and both give comparable results. 

To obtain absolute rest-frame I- and R-band magni- 
tudes for the cluster members we determine the appropri- 
ate K-correction K\^(z) for the cluster (deflector) redshift 
Zd = 0.451 elliptical galaxies and FORS1 filters using the 
GISSEL library (Bruzual A. & Chariot 1993) and obtain 
Ki{z d ) = 0.378, K R (z d ) = 0.747. In addition we apply 
galactic extinction Ai^ to the measured isophotal mag- 
nitudes and assume zero evolutionary correction. We use 
Ai = 0.121, and A R = 0.166 from NED, where the values 
are obtained from Schlegel et al. (1998) and Cardelli et al. 
(1989). 

2.2. Background galaxy catalog 

In contrast to the procedure we describe above, we use 
aperture magnitudes for the redshift determination of the 
background sources. The reason for using aperture instead 
of isophotal magnitudes is that for faint, noisy sources an 
estimate for the true object isophote is hard to achieve 
and can bias our results for these sources. The diameter 
of the aperture is set to twice the value of the seeing given 
in Table 1. In principle, one should degrade all the im- 
ages to match the seeing of the worst one (in our case U) . 
However, the effect is negligible compared to the photo- 
metric errors in the U-band, and therefore we compensate 
for that by choosing different sizes of the aperture. 

For the weak lensing analysis the R- and I-band expo- 
sures were used. As outlined in Sect. 2, the I-band serves 
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Fig. 2. The photometric redshift distributions of the cluster 
members selected as described in Sect. 2.1. 



as our primary weak lcnsing science frame. It is the im- 
age with the highest number-density of sources that can 
be used for weak lensing. Below, we cross-check our re- 
sults obtained in this band with a parallel analysis in 
the R-band. We correct all galaxies in the field for the 
PSF anisotropy and PSF smearing as described in Erben 
et al. (2001). The procedure is based on the KSB algorithm 
(Kaiser et al. 1995), in particular we use the IMCAT imple- 
mentation (http://www.ifa.hawaii.edu/~kaiser). We 
select stars from the half-light-radius vs. magnitude dia- 
gram and fit a second-order polynomial to their measured 
ellipticities. In Fig. 3 we plot the measured PSF variation 
for the I- and R-band data. 

For the final weak lensing catalogue only sources hav- 
ing photometric redshift estimate z p hot > 0.55 are con- 
sidered. We end up with N g — 210 background sources 
for the I-band data (giving 15 galaxies per arcmin 2 ), and 
with N g = 140 (10 arcmin -2 ) for the R-band. The result- 
ing redshift distributions for both catalogues are given in 
Fig. 4, the mean photometric redshift of the samples are 
(zi) = 1.18 and (z n ) = 1.14. 

3. Searching for multiply imaged candidates of 
RX J1347.5-1145 

Thus far, five arc candidates for this cluster have been re- 
ported in the literature. The first two were discovered by 
Schindler et al. (1995), and shallow HST STIS images re- 
vealed three additional ones (Sahu et al. 1998). These five 
arcs (A1-A5 as labelled by Sahu et al. 1998, see Fig. 7) 
are not all images of the same source. As is obvious from 
Fig. 6 most of them have different colours and surface 
brightnesses. Since gravitational lensing conserves both 
they belong to at least three different sources. However, 
two of these arcs (A4 and A5) do have the same colours 
and we consider them to be images of the same source. 



Although A4 has the appearance of a very straight, edge- 
on spiral galaxy (see Fig. 7), it can still be lensed, since 
the cluster members to the south west of it can produce a 
sufficiently strong tidal field to cause such a morphology. 
We note that the arc A3 considered by Allen et al. (2002) 
to belong to this system as well has different colours (see 
Table 3). Judging by the Fig. 6 one would think that Al 
and A3 are multiple images of a single source as well, 
however the figure is a composite of 3 bands only and the 
detailed photometry shows that this is not the case (see 
Table. 3). 

We detect new arc candidates using the I-band and 
Ks-band image, as well as the combined (following the 
procedure described in Szalay et al. 1999) UBVRIKs im- 
age. In individual bands some of the arcs could not be 
significantly detected. In particular, we report here on the 
discovery of a red double-component arc candidate to the 
south-west of A4, which we designate with labels Bl and 
B2. The two components formed in the middle of a con- 
centration of cluster members. Their extreme red colour 
suggests that it is either a highly reddened galaxy at z ~ 1 
or it is a galaxy at z > 5. In addition we detect in the 
vicinity of the system B a long thin arc candidate (C), 
which was also presented in Lenzen et al. (2004) as num- 
ber 3 (see Fig. 8). In the vicinity of A2 we detect addi- 
tional four arc candidates and denote them as D1-D4 (see 
Fig. 7). However, we do not claim that these components 
components belong to the same multiply imaged system, 
although their configuration is suggestive for that. Since 
these candidates are very faint, no reliable photometry can 
be obtained; the same is true for the arc candidate E. We 
use SExtractor to measure the ellipticities of these arcs 
from the I-band (systems A, C), Ks-band image (system 
B, due to its extreme red colour), and combined UBVRIKs 
image (systems D, E; since they can not be significantly 
detected in individual bands) see Table 2. We detect 
more possible arc candidates (labelled only with arrows in 
Fig. 7). They are at the limit of the detection level and 
therefore their associated errors are too large for them to 
be used for our analysis. 

Starting from the most plausible candidate multiple 
image system A4-A5 we search for additional images be- 
longing this system in an automated fashion. The aperture 
magnitudes of an image in either Nf = 6 or Nf = 5 filters 
rriij are compared with the magnitudes rrijj of all other 
images in the field (where i is in our case the index of A4 
or A5). We use the x 2 approach 



•y 2 - 



£ 



J iJ 



(2) 



where fiij is the relative magnification between the images 
i and j, and mj and ajj are the magnitude measure- 
ment errors. Since lensing is achromatic we can evaluate 
fiij by forcing Oxfj/dfiij = to hold. The resulting xfj 
function follows a ^-distribution with Nf — 1 degrees of 
freedom. The best fitting images are then further visually 



M. Bradac et al.: The cluster mass distribution of the most X-ray luminous cluster RX J1347.5— 1145 



I 



2500 

2000 
1500 

1000 E- 

500 




7t — i — I — i — r 



_l I I I I I I I I I I L 



1000 2000 
x 



R 



2500 
2000 
1500 
1000 F- 
500 E- 




f \ A , ^ 



J L 



1000 2000 
x 



0.02 



0.01 - 



o 
o 

CM 
D 



-0.01 



-0.02 



O 

u 

CM 



-0.02 





e1 corr 



0.02 



0.02 
0.01 


-0.01 F- 
-0.02 



n I r 



* x * 

**x* X 
x£x* 



-0.02 





e1 corr 



0.02 



Fig. 3. The upper panels show the spatial variation of the PSF anisotropy in I (34 stars) and R (32 stars). The length of the 
sticks give the amplitude of the stellar ellipticities (the length corresponding to |e| = 0.01 is given at the top-left). The maximum 
cllipticity is around 2%. The lower panels show the ellipticity distribution after correction with a second-order polynomial. The 
formal residuals are about 0.005 in each component. Because of the small field of view and the small number of stars, the true 
errors are probably higher but difficult to estimate. 



analysed and tested for the conservation of surface bright- 
ness. A possible counter image candidate to A4 and A5 
was found, which we label with AC. All three are encir- 
cled dashed- yellow in Fig 6, their photometric properties 
(and the properties of A1-A3) are listed in Table 3. 

Using six flux measurements in UBVRIKs for the red- 
shift determination of A4 and A5 and five in UBVRI for 
the AC (it is located at the edge of the Ks-band image 
and therefore the Ks photometry is not reliable) we find 
that A4 and A5 are consistent with being at a source red- 
shift of z s ~ 1.76. Unlike for other background objects (see 
Sect. 2.2), we use isophotal magnitudes to obtain reliable 
rcdshifts for A4 and A5 here due to the large ellipticity of 
the arcs. The redshift determination is in agreement with 
Ravindranath & Ho (2002) who, based on the absence of 
the 0[II] line in their spectrum, predict the redshift of 
A4 to be > 1.04. The redshift estimate of AC using 5 fil- 
ters is 1.3; however, also the redshift estimates of arcs A4 
and A5 are 1.3 if we use only 5 filters. All three prob- 
ability distributions for the redshift estimates are very 
broad and the higher redshift of 1.76 is consistent with 
the photometric data in all three cases. In the redshift 
regime 1.2 < z < 2 the main features in the spectral en- 



ergy distribution (Lyman break, Balmer break, etc.) lie 
outside of the optical bands and therefore the NIR pho- 
tometry is important. We therefore use the estimated pho- 
tometric redshift from UBVRIKs of z s ~ 1.76 from now 
on. Unfortunately, the redshift estimate for the multiply 
imaged system can substantially influence the combined 
cluster mass reconstruction (the position of the critical 
curve changes with redshift). We investigate this effect in 
Sect. 4.3. 



Within the errors, the three images have the same 
colours as well as the same surface brightnesses (see 
Table 3). In addition, the photometric redshift estimate 
(using 6 filters) is the same for A4 and A5. The colours 
and peak surface brightnesses of the counter image are 
also consistent, however due to its smaller apparent size 
its photometry is less reliable. There are more candidate 
multiple image systems in this field; they will be the sub- 
ject of a future study. 



M. Bradac et al.: The cluster mass distribution of the most X-ray luminous cluster RX J1347.5— 1145 



o 
oo 



o 

CD 



r ° 



O 



o 



o 



'phot 



(a) 



'phot 



(b) 



Fig. 4. The redshift distributions of background sources used for weak lensing analysis (only sources with z phot > 0.55 were 
considered) for I-band a) and R-band catalogue b). The mean photometric redshift of the samples is (zi) — 1.18 and (zr) = 1.14. 



4. Cluster mass reconstruction of 
RX J1347.5-1145 

In this section we present the mass modelling of the cluster 
RX J1347.5-1145. We first give a short outline of the 
method, a full account of it can be found in Paper I. 

Table 2. The properties of the arcs A1-A5 and the candidate 
counter-image AC used in the strong lensing analysis of the 
cluster. We also present additional arc candidates (systems B, 
C, D, and E) - see also Fig. 7. The properties of systems A, 
and C are measured from the I-band, while B is measured 
from the Ks-band image, and D and E are measured from the 
combined UBVRIKs image. The positions and position angles 
are given with respect to the brightest cluster member; the 
position angle of 90° means that the arc is tangentially aligned 
with the BCG. All are measured with SExtractor. Because of 
the proximity of Bl and B2 to a cluster member we can not 
measure their ellipticities accurately. 



Arc 


8i [arcmin] 


82 [arcmin] 


N 


PA[deg] 


Al 


-0.3460 


0.4547 


0.571 


82.0 


A2 


-0.1253 


0.5103 


0.505 


79.1 


A3 


0.6644 


0.0090 


0.613 


96.1 


A4 


0.3314 


-0.4980 


0.713 


99.8 


A5 


-0.2891 


-0.7009 


0.333 


99.4 


AC 


-1.0241 


0.6509 


0.327 


65.8 


Bl 


0.7440 


-0.6873 






B2 


0.6971 


-0.7653 






C 


0.3299 


-0.4985 


0.735 


99.8 


Dl 


0.1985 


0.4492 


0.219 


81.8 


D2 


0.3023 


0.3688 


0.283 


80.31 


D3 


0.3812 


0.2562 


0.364 


56.5 


D4 


0.4771 


0.1162 


0.488 


106.5 


E 


-0.4048 


-0.3581 


0.552 


108.97 



4.1. Short outline of the method 

The main idea behind the method is to parametrise the 
cluster mass-distribution by a set of model parameters, 
where this parametrisation is chosen as generic as possi- 
ble. In our case we use the gravitational potential ip on a 
regular grid. We factorise the redshift dependence by the 
so-called "cosmological weight" function Z(z), as defined 
in Paper I (see also Bartelmann & Schneider 2001). 1 
We define the x 2 -function 

x 2 (#) = Xttyk) + vR&k) + xliWk) ■ (3) 

where Xe(V'fe) is the contribution from weak lensing and 
XmCV'/c) from strong lensing. In addition, the regularisa- 
tion R(ipk) with regularisation parameter rj is employed in 
order to penalise any models that would follow the noise 
pattern in the data. We minimise the \ 2 function with 
respect to ipk by solving the equation dx 2 (ipk)/dijjk = 0. 
This is in general a non- linear set of equations, and we 
solve it in an iterative manner. We linearise this sys- 
tem and starting from some trial solution (i.e. k^°\ 
and a(°>) we repeat the procedure until convergence is 
achieved. We showed in Paper I that different models used 
as a trial solution do not influence results significantly. We 
try to confirm this result by investigating different models 
here as well. 

4.2. Initial conditions for the method 

For the purpose of obtaining the initial values for 

and we first investigate the signal from the aver- 
aged tangential ellipticities and fit these using the singular 

1 To evaluate the angular diameter distances we assume the 
ACDM cosmology with Q m = 0.3, Qa = 0.7, and Hubble con- 
stant H = 70km s _1 Mpc -1 . 
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isothermal sphere SIS model (hereafter called IS* scenario). 
The tangential ellipticities are given by 

e t = -3?[ee- 2i ^] (4) 

where <j) specifies the direction to the source galaxy with 
respect to the the brightest cluster galaxy (BCG). In Fig. 5 
we plot the average tangential ellipticities versus projected 
radius (9) in radial bins centred on the BCG, containing 
50 (35) galaxies each. Both, the I- and R-band data give 
comparable results. We note that the tangential ellipticity 
signal is still high at the edge of the field (e t ) ~ 0.1, thus 
making this data inaccessible for standard weak lensing 
techniques aiming to determine the mass, since on this 
relatively small field one cannot break the mass-sheet de- 
generacy by simply assuming k ~ at the field edges. 

We fit an SIS profile to the individual tangential ellip- 
ticities (not binned), the model ellipticities are calculated 
using the redshifts of these sources. The resulting line-of- 
sight velocity dispersion is cri sis = 950 ±60 kms -1 for the 
I-band data and 0r ; sis = 960 ± 70 kms -1 for the R-band 
(both ltr error bars). The tangential ellipticity as a func- 
tion of radius for this model is plotted in Fig. 5 (dashed 
line) for the average source redshifts of (z\) = 1.18 and 
(zr) = 1.14 (see Sect. 2.2). In addition, the absence of the 
lens is excluded with more than lOc significance in both 
bands (all minimisations and error analysis in this subsec- 
tion are performed using C-mimiit from James & Roos 
1975). 

The line-of-sight velocity dispersion estimates are 
higher than the measured velocity dispersion from Cohen 
& Kncib (2002), and lower than previous weak lensing, 
strong lensing and X-ray measurements. However, in the 
optical it is evident that the cluster has a lot of structure 
and therefore the SIS profile does not describe the cluster 
adequately. It has at least two main components; in ad- 
dition there is X-ray emission off-centred from the BCG. 
Furthermore, at the scales where we measure the profile, 
< 400/i -1 kpc, the profile of the cluster is probably not 
isothermal (see e.g. Navarro et al. 2004). Therefore, the 
values of a obtained in this manner should not be trusted, 
we only use them for one of the initial models for 

Another possibility to obtain initial conditions is to use 
the multiple image information for the cluster. We perform 
a very rough analysis by using the data for the arc system 
A4-A5-AC (given in Table 2). In addition to the image 
positions we also use image ellipticities as constraints. The 
model consists of a non-singular isothermal ellipse (NIE) 
(Keeton & Kochanek 1998), given by 

k{6') = bo , (5) 

2 V^( rc2nis + W)2 ) + ( ^ )2 

where 6' is calculated w.r.t. the major axis of the cluster 
surface mass density. We allow the centre of the cluster 
potential c i, the scaling bo, ellipticity |e g |, and the posi- 
tion angle <p g to vary. Following the prescription of Kncib 
et al. (1996) we also include the 10 brightest cluster mem- 
bers in the I-band to the mass model. They are modelled 



as non-singular isothermal spheres with a line-of-sight ve- 
locity dispersion er n j s and core radius r c , n i s following 

cr n is oc L 1/4 , r C! n is oc L 1/2 , (6) 

The proportionality constants were chosen such that the 
I-band magnitude mj = 17.5 galaxy would have cr sis = 
300 kms" 1 and r c , nis = O'.l (the BCG has mj = 17.8). 
We also fix the core radius of the cluster to 0'.3. These 
constants are not allowed to vary. The best fit model 
for this system has values of {6 c i,i, &c\.2, bo, |e g | , cf> g } = 
{-0'.21,-0'.10,0'.97,0.3,0.8}. 

We stress here that it was not our aim to obtain a de- 
tailed strong lensing cluster-mass model, since it will only 
be used for the initial values of reconstruction. The mul- 
tiple image system used here is independently included 
in the non-parametric reconstruction. We have shown in 
Paper I (and also confirm this in Sect. 4.3) that the re- 
construction depends little upon the details of the initial 
model we use; for this reason a detailed modelling is not 
needed. In particular, the precise choice of those param- 
eters that we did not vary in the modelling is not very 
relevant in our case. For the same reason we also do not 
include additional multiply imaged candidates in the anal- 
ysis. 




Fig. 8. A ~ l'.O x 0'.6 cutout of the Ks-band image of 
RX J1347.5— 1145 showing one of the arcs used for strong lens- 
ing, the newly discovered very red arc system B1-B2 (two im- 
ages almost merging, c.f. Fig. 6) and the long thin arc C. 



4.3. Combined weak and strong lensing mass 
reconstruction of RX J 1347. 5- 1145 

We apply the mass-reconstruction method to the I- and 
R-band data of RX J1347.5— 1145 and the strong lensing 
system A4-A5-AC only (see Table 2). We use three differ- 
ent initial models for ; IM is the best fit model from the 
strong lensing analysis of the cluster, the IS model is the 
best fit SIS model to binned tangential ellipticities (cen- 
tred on the brightest cluster member) - both presented in 
Sect. 4.2 - and 10 has «W = fj = 0, q<°) = 0). The 
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Fig. 5. Average tangential ellipticity (et) vs. projected radius {8} in radial bins centred on the brightest cluster member contain- 
ing 50 galaxies per bin for the I-band data a) and 35 galaxies for the R-band data b). The errors are obtained by randomising 
the phases of the measured ellipticities, while preserving their absolute values. 100 randomisations were performed. The dashed 
line shows the best-fit SIS profile to the unbinned data, for the I-band we obtain d.sis = 950 ± 60 kms -1 and for the R-band 
or, sis = 960 ± 70 kms -1 (both la error bars). It is plotted here for the average source redshifts of (zi) = 1.18 and (zn) = 1.14. 
In the interval 0' < (8) < 0'.25 the expectation value of the observed ellipticity (for these models) is given by l/g(0, z)* and the 
tangential ellipticity profile has a steep gradient at {8) ~ O'.l (for both bands) where Z({zr.,i))k = 1 [since g — Zy/ (1 — Zn)]. 
At (8) > 0'.25 the expectation value of the observed ellipticity is given by g(6, z). 



initial regularisation parameter is set to r\ = 400 for the I- 
band and rj = 200 for the R-band. It is adaptively adjusted 
in each iteration step such that the resulting X 2 /A^ g ~ 1. 
The resulting K-maps are given in Fig. 9. We also overlay 
the contours from Fig. 9al to the colour composite image 
in Fig. 6. 

We estimate the mass within the cylinder of a ra- 
dius of 1'.5 (for the cluster redshift = 0.451 this corre- 
sponds to 360/i _1 kpc), the estimates are given in Table 4. 
The projected mass of the cluster is estimated to be 
M(< 360 fc^kpc) = (1.2 ± 0.3) x 10 15 M Q . The error 
was estimated by bootstrap resampling the background 
galaxies in the weak lensing catalogues. This means that 
for each catalogue we randomly select N g galaxies with 
replacement, if a galaxy is selected twice (or more) we as- 
sign double (or multiple) weight to its Xe contribution. 
We generate 10 new catalogues and perform a new mass 
reconstruction; the error is then given as the variance of 
these estimates. It is larger than what we obtain from sim- 
ulations in Paper I, which is partly attributed to the fact 
that we only use a three-image and not a four-image sys- 
tem here. However, within the given errors the results for 
both bands and for different initial models are consistent. 

The projected mass from XMM measurements (Gitti & 
Schindlcr 2004) within a cylinder of the same radius as we 
use is given by M x (< 360ft _1 kpc) = (0.7 ±0.2) x 1O 15 M 
(Myriam Gitti, private communication). The resulting 
mass from the strong and weak lensing mass reconstruc- 
tion is higher and marginally consistent with X-ray mea- 
surements. If the mass estimate is extrapolated at larger 



radii (assuming an isothermal profile) it is also consistent 
with the previous weak-lensing results by Fischer & Tyson 
(1997). It is however significantly larger than the mass es- 
timate obtained by the velocity dispersion measurement 
of Cohen & Kneib (2002). 

A possible explanation for the discrepant dynamical 
mass estimates was presented by Cohen & Kneib (2002). 
They argue that the cluster is most likely in a pre-merging 
process (with clumps merging preferentially in the plane of 
the sky). In such a scenario, until the merging is complete 
and the cluster is virialised, the dynamical cluster mass 
will be largely underestimated. On the other hand the 
X-ray temperature can be increased in such merging pro- 
cesses (thus the mass would be overestimated) and for this 
reason the south-east quadrant is excluded (the surface 
brightness profile is determined by averaging data only 
in the other three quadrants) in the X-ray analysis. The 
temperature measurements from Gitti & Schindler (2004) 
thus further supports the merger hypothesis. However if 
there is some extra mass present in the excluded quadrant 
(as suggested by our mass maps), the mass estimate ob- 
tained in this way from X-rays will be underestimated. If 
the hypothesis is correct, traditional mass estimates rely- 
ing on equilibrium assumptions fail and gravitational lens- 
ing (with high quality data) provides the most accurate 
estimate for the cluster mass. 

We note, however, that our results depend upon the 
correct redshift determination and identification of the 
members of the multiple image system we use. If we put 
the multiple image system to a redshift of ~ 3 (~ 1.3), 
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Fig. 6. The BRK colour composite of the ~ 3 x 3.2 amin 2 field of RX J1347.5-1145. Overlaid is in white contours the combined 
weak and strong lensing mass reconstruction from Fig. 9al. The contour levels are the same (the field here is smaller), we 
smooth them here using a Gaussian kernel characterised by a — 5" for clarity of the plot. Yellow circles show the multiple image 
system we use. North is up and East is left. 



the estimated mass decreases (increases) by ~ 10%. If the 
images do not belong to the same system, the changes 
might be even more drastic. However, at least for the two 
arcs A4 and A5, based on their photometric properties, 
we consider this possibility less likely (see Sect. 3). As a 



test we have also performed the reconstruction using only 
the two arcs A4 and A5. The results remain unchanged, 
however the scatter between the three initial models and 
the errors are larger by a factor ~ 2. 
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Fig. 7. The 2.7 x 2.4 arcmin 2 combined UBVRIKs image of the RX J1347.5— 1145 showing previously known (system A and C) 
and newly discovered (systems B, D, and E) arc candidates (see Table 2). All are marked with a label and an arrow, further arc 
candidates are marked with an arrow only. The image was scaled non-linearly and the scaling varies across the field (in order 
to enhance the details of the image). 



Further, the results rely on the correct determination 
of the photometric redshifts for the weak lensing sources. 
The random error of the determination is not crucial, the 
problem are the systematic uncertainties. It is not ex- 
cluded that e.g. some foreground sources get assigned a 
high redshift and thus diluting (if they are randomly ori- 
ented) or enhancing the signal (if they are aligned). In 
addition, outliers can have Z(z) assigned which is very 
different to their real cosmological weight. These outliers 
were considered in Paper I, they were chosen at random 
and their fraction was taken to be 10%. Still, their pres- 
ence did not significantly change our conclusions. If, how- 
ever, their fraction is higher and/or more importantly if 
they bias the final redshift distribution, this can bias our 
mass estimate. 

An additional test for the accuracy and reliability of 
our model could be performed by using its predictive 
power. Namely, if the model is well constrained it should 
be capable of predicting the position of e.g. the counter 
image to the arc Al (providing its redshift determination 



is correct). We have tested our models using the follow- 
ing procedure. Using the resulting potential from strong 
and weak lensing reconstruction we project (using bilin- 
ear interpolation and finite differencing) the position 9 S of 
an arc candidate (e.g. Al) back to the source plane and 
denote the resulting position as y s (0 s ). Then we search 
for all possible solutions 9 satisfying the non-linear set of 
equations y(9) — y s (9 s ) = 0. These should then lie close 
to the possible counter image candidate(s). However, since 
our model is tightly constrained only in the vicinity of 
multiple images we use (A4, A5, and AC), the scatter of 
possible solutions is large. This issue could however be eas- 
ily resolved in the future with e.g. ACS observations, since 
many more arc candidates will be found and their mor- 
phology can be obtained allowing for unambigous identifi- 
cation of multiple imaged systems and tighter constraints 
of the model. 
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Fig. 9. K-maps obtained from combined strong and weak lensing reconstruction of the cluster RX J1347.5— 1145. Left panels 
show the reconstructions using I-band data and for the ones on the right we use the R-band data. For the I-band data we have 
N g — 210 background source galaxies, and for the R-band N g = 148, all with known photometric redshifts. In al) and a2) 
we use the best fit model from the strong lensing analysis of the cluster IM as initial condition, in bl) and b2) we use the IS 
model (SIS model fitted to tangential ellipticities, centred on the brightest cluster member) - for both IM and IS see Sect. 4.2. 
In cl and c2) we use 10, i.e. k^ ' = on all grid points. The positions of the 10 brightest cluster members are plotted as white 
circles. 



4.4. Rest-frame I- and R-band brightness distribution 
and mass-to-light ratio of RX J1347.5-1145 

To obtain the cluster brightness distribution and aperture 
luminosity we proceed as follows. Using colour and redshift 



information for the selection criteria described in Sect. 2.1 
we determine the luminosities of the cluster members in 
the field. They are smoothed using a Gaussian kernel char- 
acterised by (7 — 9", resulting in the brightness distribu- 
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Fig. 10. The I-band a) and R-band b) brightness distribution of the RX J1347.5-1145 in 1O 12 L / arcmin 2 . The cluster 
members are selected using the colour cuts described in the text, the ten brightest ones (in I) are plotted as white circles. Their 
luminosities of the cluster members have been smoothed using a Gaussian kernel characterised by a — 9". 



tion shown (for I-band only) in Fig. 10. We then determine 
the aperture luminosity L a by adding the luminosities of 
the cluster members within 360/i _1 kpc radius around the 
BCG. 

The resulting I- and R-band aperture luminosities 
are L a ,i(< 360/i _1 kpc) = 3.1 x lO 12 i and L a>R (< 
360ft- _1 kpc) = 2.2 x 1O 12 L , respectively. The mass-to- 
light ratios (M/L) are M/L Y = 400 ± 15OM /Li i0 and 
M/L R = 550 ± 15OM /L R , Q . The cluster has only 300 
members across the observed field, it is under-luminous in 
optical bands. In addition, we are measuring the M/L ra- 
tio in the inner part of the cluster, which might not reflect 
the M/L ratios measured out to ~ 1 Mpc distances from 
cluster centres usually quoted in the literature. 

The first concern with luminosity estimates is com- 
pleteness. For this purpose we fit the Schechter lumi- 
nosity function (Schechter 1976) to the cluster member 

Table 4. Reconstructed mass of RX J1347.5— 1145 within a 
cylinder of 360fe _1 kpc radius around the BCG from I-band 
(left) and R-band (right) weak lensing data and one candidate 
3-image system. Three different models have been used. 
We use the best fit model from the multiple image system IM, 
IS is the best fit SIS model from the process of parametrised 
fitting of weak lensing data and 10 has = on all grid 
points (see Sect. 4.2). In brackets we give for comparison the 
velocity dispersion of an SIS having the same enclosed mass 
within 360/i _1 kpc. 





Mi 


[cti.sis] 


Mr 


[or, sis] 




[1O 15 M ] 


[kms- 1 ] 


[1O 15 M ] 


[kms" 1 ] 


IM 


1.37 ±0.04 


[1900] 


1.31 ±0.03 


[1860] 


IS 


1.2 ±0.1 


[1800] 


1.1 ±0.1 


[1700] 


10 


1.1 ±0.1 


[1700] 


1.1 ±0.1 


[1700] 



counts as a function of absolute magnitudes M\ and Mr.. 
The resulting best-fit characteristic magnitudes are M{ = 
-22.4±0.2 and M K = -21.7±0.3 and the faint end slopes 
are given by cm = -0.89 ± 0.06, and a R = -0.9 ± 0.1. We 
conclude that our catalogues are complete to a magnitude 
Mj* R — 4 and therefore the contribution from incomplete- 
ness is negligible. 



A more severe concern is the contamination by non- 
cluster members and rejection of the actual members. In 
order to check against this, one needs to investigate the 
galaxy population "outside" the cluster region (on images 
taken with the same photometric conditions and depth 
as the images we use). A slightly different approach for 
the purpose of the M/L calculations can be followed by 
defining an outer aperture at the edge of the image and 
subtracting the luminosity density in that aperture from 
that in the inner portions of the cluster. The same ap- 
proach needs to be undertaken when calculating the mass. 
If the M/L is constant across the field, this would give 
its correct value. Unfortunately our observed fields span 
only ~ 450 kpc h^ 1 around the brightest cluster galaxy 
and therefore this approach is not reliable. We conclude 
that the error budget on luminosity is dominated by the 
systematics of the cluster member selection and contami- 
nation and is very difficult to estimate. We investigate two 
different selection criteria for cluster members in Sect. 2.1; 
we used colour information as well as the photometric red- 
shifts. The aperture luminosities from these two criteria 
are consistent at the 5%-level. These two approaches share 
similar systematics, both use the same magnitude mea- 
surements, and for the colour-colour selection blue galax- 
ies are added using photometric rcdshift measurements. 
However, in order to estimate the M/L the mass determi- 
nation is a dominant source of error. 
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5. Conclusions 

The case of RX J1347.5— 1145 has been a cause of many 
puzzles in the past. Very discrepant mass estimates are 
given in the literature, and unfortunately this cluster is not 
the only case where the mass measurements have proven 
to be difficult. We have applied a new mass reconstruc- 
tion method to deep optical data using a multiple-image 
system with three images selected based on their colours 
and redshifts. Our main conclusions are the following. 

1. The combined strong and weak lensing mass recon- 
struction confirms that the most X-ray luminous clus- 
ter is indeed very massive. If the rcdshift and identi- 
fication of the multiple-image system, as well as red- 
shifts used in weak lensing data, are correct we esti- 
mate the enclosed cluster mass within 360 /i _1 kpc to 
M{< 360/1"^) = (1.2 ±0.3) x 1O 15 M . 

2. The reconstruction shows a south-east mass exten- 
sion, compatible with the X-ray measurements (see e.g. 
Gitti & Schindlcr 2004, Allen et al. 2002). 

3. A single SIS fit to the average tangential ellipticities 
does not give a reliable estimate for the enclosed mass 
within 360/i _1 kpc; detailed modelling needs to be per- 
formed. 

4. We have demonstrated the feasibility of breaking the 
mass-sheet degeneracy in practice by using shape mea- 
surements and adding the information on individual 
redshifts, without any assumptions regarding the clus- 
ter potential. 

In addition we measured the corresponding mass-to-light 
ratio of the cluster within 360 /i~ 1 kpc. We find that the 
cluster is more luminous in the rest-frame I-band than R- 
band, which is expected due to the presence of many old 
(red) elliptical galaxies in clusters. The resulting mass-to- 
light ratios are high, both in rest-frame I- and R-band, 
giving M/Li = 400 ± 15OM /Li i0 and M/L R = 550 ± 
15OM q /Lr,,0. These values are higher than typical values 
for clusters claimed in the literature (~ 200 for R-band). 
However it is difficult to compare our results with existing 
measurements of the mass-to-light ratios, since they are 
usually performed at larger radii not accessible with our 
data. 

In the course of this research wc discovered one new 
extremely red arc candidate (system B) at ~ 1 arcmin dis- 
tance from the BCG. Unfortunately its redshift can not 
be measured, as it is significantly detected only in the Ks 
band. Further arc candidates are discovered from the com- 
bined colour image, suggesting that the cluster is indeed 
very centrally concentrated. In addition, the enclosed mass 
we obtain using the combined reconstruction also fits rea- 
sonably well the standard mass vs. X-ray luminosity rela- 
tion (see Rciprich & Bohringer 2002), provided we assume 
the model to be isothermal (which for the same enclosed 
mass as our reconstructed model means a ~ 1800 kms -1 ) 
to a radius of r 2 oo frequently used to determine the rela- 
tion. 



The mass-reconstruction of RX J1347.5— 1145 can be 
significantly improved. Deep HST imaging would greatly 
help in identifying and confirming new multiple-image sys- 
tems, thus allowing more detailed modelling. In addition, 
not only the centre of the light for each of the arcs can be 
used as constraints, but also their morphology. As men- 
tioned in Paper I, the reconstruction technique with adap- 
tive grid at image positions can be used for these purposes. 
Further, spectroscopic redshifts need to be obtained for 
the multiple-image system candidates as well as for the 
cluster members (to obtain velocity dispersion measure- 
ments from a large sample). Deep, wide- field imaging data 
of this cluster will help us to improve the weak lensing con- 
straints also at larger radii than presented here. A large 
number density of sources that can be used for weak lens- 
ing accessible by ACS (> 120 arcmin^ 2 ) would greatly 
improve the accuracy of the mass estimate and enable us 
to resolve substructures in the cluster. The details of the 
reconstruction can be used to reliably determine the clus- 
ter profile. 

In conclusion, even without the best data quality that 
can be reached at present, we were able to perform a de- 
tailed cluster-mass reconstruction of the most X-ray lu- 
minous cluster RX J1347.5— 1145. The method has also 
shown a high potential for the future. If the highest qual- 
ity data is used, a combination of strong and weak lensing 
has proven to offer a unique tool to pin down the masses 
of galaxy-clusters as well as their profiles and accurately 
test predictions within the CDM framework. 
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Table 3. The photometric properties of the arcs Al, A2, A4, A5, and the candidate counter image AC. Given are three colours 
(tub ~ mi, my — mi, and rriR — mi) in magnitudes (measured from the isophotal magnitudes), VRI peak surface brightnesses 
SV,r,i (in magnitudes), and photometric redshifts. For A1-A5 we determine them using 6 bands, for AC Ks is not available. If 
objects belong to the same source the colours and surface brightnesses need to be conserved. 





ma — mi 


mv — mi 


ma — mi 


Sv 


Sn 


Si 




Al 


2.20 ±0.07 


1.78 ±0.05 


0.98 ±0.05 


24.37 


23.75 


22.72 


0.69 


A2 


3.66 ± 0.09 


2.51 ±0.05 


1.38 ±0.05 


24.16 


23.44 


21.88 


0.73 


A3 


1.52 ±0.07 


0.96 ±0.05 


0.54 ±0.05 


24.23 


23.89 


23.19 


1.65 


A4 


0.99 ± 0.07 


0.81 ±0.05 


0.53 ±0.05 


23.60 


23.30 


22.60 


1.76 


A5 


1.08 ± 0.07 


0.88 ±0.05 


0.57 ±0.05 


23.60 


23.29 


22.60 


1.70 


AC 


1.28 ± 0.07 


0.88 ±0.05 


0.52 ±0.05 


23.90 


23.50 


22.98 


1.30 a 



a The redshift of AC was determined using 5 bands only. It is consistent with redshifts of A4 and A5 if they are also determined 
without Ks-band information. 



